Scaling of the entanglement spectrum near quantum phase transitions 
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The entanglement spectrum describing quantum correlations in many-body systems has been 
recently recognized as a key tool to characterize difTerent quantum phases, including topological 
ones. Here we derive its analytically scaling properties in the vicinity of some integrable quantum 
phase transitions and extend our studies also to non integrable quantum phase transitions in one 
dimensional spin models numerically. Our analysis shows that, in all studied cases, the scaling of the 
difference between the two largest non degenerate Schmidt eigenvalues yields with good accuracy 
critical points and mass scaling exponents. 

PACS numbers: 



The frontiers of condensed matter physics have broad- 
ened substantially beyond the paradigms of the twentieth 
century, namely conventional superconductivity (BCS) 
and the band theory, with the discovery of topological 
quantum phases. These novel states are widely thought 
to have great potential for technological applications, 
fndccd, first attempts in this direction are already in 
progress. Some examples are the use of the edge states 
for topological quantum computations [l[ or the cheap 
synthesis of high-Tc superconductors via the doping of 
certain spin liquids Q. However, for the sake of appli- 
cations a deeper theoretical comprehension is required. 
Yet another milestone has been the gradual convergence 
of condensed matter physics and quantum information 
theory (QIT) motivated by the fact that quantum many- 
body systems arc a natural territory of quantum entan- 
glement, which is the basic computational resource in 
QIT. Such theoretical developments run in parallel and 
are stimulated by the spectacular experimental progress 
achieved in the areas of ultracold physics and condensed 
matter. New platforms such as ion traps, optical lattices, 
bosonic and fcrmionic atomic degenerate gases Q, new 
superconducting materials, exotic quantum Hall systems 
are intensively investigated as quantum simulators. On 
these devices a huge set of new phases is postulated to 
be synthesizable. In the light of these new exciting de- 
velopments, the problem of how to characterize, detect 
and simulate new quantum phases of matter remains as 
one of the greatest challenges in the field. 
Traditionally, the characterization of quantum phases 
and their phase transitions has been based on local order 
parameters Q, and the response to linear perturbations 
given by low order correlators e.g. (SiSj), where Si is 
an observable of site i in a lattice system. Such descrip- 
tion is in accordance to the "standard" Ginzburg-Landau 
scenario of phase transitions. Here the order, which is as- 
sociated to the breaking of some symmetry, is manifested 
by a change in the expectation value of a local order pa- 
rameter Q. However, there are important cases where an 
order parameter is not available, as in presence of topo- 
logical phases or more simply whenever a local parameter 



is intrinsically difficult to construct or measure. Notable 
examples arise for instance in the presence of deconfined 
criticality, like at the transition semimetal-insulator on 
the graphene or in the Haldane-Shastry model [1, or 
in long range interacting systems @. A similarly chal- 
lenging situation realizes whenever the phase diagram of 
the material under investigation is not known at all, as 
for some high-Tc superconductors or for various lattice 
tight-binding models (see for instance 0)- In all these 
cases a major problem is to identify a quantity able to 
detect a certain phase or phase transition. 
In recent years, concepts from QIT have helped to un- 
derstand the structure of quantum many-body systems 
and the computational power needed to simulate them 
0. A major step in this direction was the study on 
how the bipartite entanglement of a block A of a many- 
body system, as measured by the entanglement entropy 
S = — (0/1 log PA {pA being the reduced density matrix), 
scales with the size ha of the block Later on, it 
was rigorously proven that for all ID gapped quan- 
tum systems described by short range Hamiltonians, the 
entanglement entropy saturates to a constant indepen- 
dently of the size of the block. This behavior of the en- 
tropy, encountered in several areas of physics, is termed 
"area law" At the quantum phase transition (QPT) 
instead the entanglement entropy diverges logarithmi- 
cally with the block size as S* ~ c\og{nA) [13 - [l4l | where c 
is the central charge of the conformal field theory (CFT) 
describing the critical point. Out of criticality, the many- 
body system is not any longer entirely constricted by the 
central charge but has a dependence on the specific per- 
turbations tuned to move from the transition point. It is 
in this regime that one expects the entanglement spec- 
trum, i.e. the set of the eigenvalues {A^} of the reduced 
density matrix pA, to contain information that is not in- 
cluded in the entanglement entropy, being this quantity 
a single number. Indeed Li and Haldane [isj and im- 
mediately later various authors [l^ pointed out that the 
entanglement spectrum of topological systems, after a cut 
in a suitable Hilbert space, shows peculiar features in de- 
generacy, gap, distribution and it can be a valuable tool 
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for the investigation of many-body systems. The struc- 
ture of the entanglement spectrum reflects some specific 
features of the phase in which the system is, like symme- 
tries, edge states, etc. In view of these developments, we 
have recently examined 17[ the entanglement spectrum 
focusing on its specific scaling features near some QPTs. 
In particular, using finite size scaling (FSS) we have 
numerically demonstrated that in the longitudinal and 
transversal Ising model, the Schmidt gap, i.e., the differ- 
ence between the two largest non trivially degenerated 
eigenvalues of the reduced density matrix A A = Ai — A2, 
correctly signals the critical point and scales with criti- 
cal exponents related to the CFT describing the transi- 
tion point itself and to its specific perturbations. Recent 
studies in 2D systems also show the scaling of the entan- 
glement spectrum near phase transitions jisl [l9j . 
Here, we extend our previous results first by deriving 
analytically the scaling of the Schmidt gap for spin- 
1/2 Ising model with transverse or longitudinal magnetic 
field. Second, we analyze numerically the scaling proper- 
ties of the entanglement spectrum in non integrable cases. 
Our main result is that for both, integrable and non- 
integrable, models the Schmidt gap always shows scaling 
behavior and the FSS of this quantity yields the mass 
scaling exponent and the precise critical point. 
We stress that for a lattice system there are two differ- 
ent concepts of integrability: the first one concerns lattice 
integrability, equivalent to the total solvability of the lat- 
tice model (for details see [13, HH), while the second one 
involves the solvability of the infrared continuos theory 
describing the low energy sector of it (the prominent one 
close to criticality) [2l| . When not explicitly specified, in 
the following we will mean the second type of integrabil- 
ity. 

The paper is organized as follows. In SectionU we first re- 
view previous results concerning the entanglement spec- 
trum and the Schmidt gap at criticality for finite size 
systems. Then, we derive analytically the closing of the 
Schmidt gap in the scaling regime (out of criticality) us- 
ing the Baxter and Cardy's conjecture that for integrable 
models near criticality one can connect the entanglement 
spectrum to the characters of the conformal field the- 
ory (CFT) describing the critical point [lO, [H-il] . Fi- 
nally, we compare the derived analytical expressions with 
the results obtained by solving the spin- 1/2 Ising chain 
in the transverse and longitudinal models. Notice that 
these two phase transitions share the same critical point 
but correspond to different universality classes. This is 
clearly reflected in values of the scaling exponents ob- 
tained for them. In Section [II] we focus on more complex 
models, including also non integrable ones. Indeed inte- 
grability is known to constrain greatly the dynamics of a 
system [21j and it could be conjectured to be responsible 
of the peculiar behavior of the Schmidt gap previously de- 
scribed. However, we observe that this is not the case. To 
this aim, using the density matrix renormalization group 
(DMRG) algorithm and the FSS [Ulil, we focus in par- 
ticular on the spin-1 Heisenberg model with an uniaxial 



anisotropy D [29|, |30[ which exhibits a Haldane topolog- 
ical phase surrounded by other quantum phases (Neel, 
dimer, large D and critical) belonging to very different 
universality classes. The Haldane-Neel transition is in 
the transverse Ising universality class, then the closing 
of the Schmidt gap reproduces the results for the spin- 
1/2 transverse Ising model. We discuss then the scal- 
ing properties in Haldane-dimer QPT described by the 
Wcss-Zumino-Witten SU{2)2 conformal theory. Finally, 
we analyze the Haldane-large D transition where the crit- 
ical exponents change along the critical line (having cen- 
tral charge c = 1) and the Neel-dimer phase transitions, 
where the central charge changes continuously along the 
transition line. Using FSS techniques, for all the above 
transitions we derive the mass scaling exponent and the 
critical point. Finally, in Section IlIII we summarize our 
results and list some open problems. 



I. THE SCHMIDT GAP IN INTEGRABLE 
MODELS: THE SPIN-i ISING CHAINS 

The entanglement spectrum is the set of the eigenval- 
ues of the reduced density matrix, known as Schmidt 
eigenvalues, for any of the two blocks obtained after a 
bipartite splitting of a system in real space. For transla- 
tional invariant systems and sufficiently large blocks, the 
point at which the bipartition is done is irrelevant, for 
simplicity we assume to be in the middle of the chain [3l| . 
After the bipartite splitting, one can always Schmidt de- 
compose any pure state, e.g. the ground state, of the 
entire many-body system as 



(1) 



where Ai > are the Schmidt coefficients with respect to 
the partition sorted in decreasing order, while and 
|0f ) are the Schmidt eigenvectors corresponding to the 
right (left) subsystem. 

The starting point in deriving the entanglement spec- 
trum is the fact that the reduced density matrix pfj of 
a semichain can be expressed as the n-th power of the 
corner transfer matrix A [2^ constructed on a certain 
portion of the classical 2D lattice associated to the quan- 
tum chain: pn = A"- [l^, [1^. For a square lattice, one 
can manage the construction such to have n = 4. If 
the lattice is integrable, A can be in turn expressed as 
A ~ exp[— where if is a fictitious non local Hamilto- 
nian defined on a semi-chain [2^ [26l | . This property also 
holds for systems that are integrable only in their con- 
tinuos infrared formulation [25j . The eigenvalues of H, 
£,k, are related to the Schmidt eigenvalues Afc as follows: 
Afe = Ce""^*^, where C is an overall constant such to ful- 
fill the constraint J^k = ^- Clearly, the Hamiltonian H 
must be gapless whenever p fj is calculated for a given sys- 
tem at a critical point. Consequently, in this point {S,k} 
and {Xk} must tend to a continuous distribution having 
statistics as in ([J) below, therefore the Schmidt gap must 
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tend to zero. We believe this property to hold also for 
nonintegrable models, since they arc continuously related 
to integrable ones in the space of the Hamiltonian param- 
eters |2l| whenever they share the same critical point. 
Let us start by analyzing first the Schmidt gap at criti- 
cality for finite size systems. In the CFT framework, the 
reduced density matrix of semichain having size (. can be 
expressed as [H, [S^ : 



1 



PR 



ZR{q) 



(2) 



where c is the central charge, Lq is the zero-generator of 
the chiral Virasoro algebra, Zii{q) = Tvq^° is the par- 
tition function of the subsystem with size £ of a torus 
(equivalent to a cylinder) [lH and q = exp(j27rT) is 
the modular parameter. The relation r ~ iK/log{£/ri) 
holds, where 77 is a regularization cutoff and k is a pos- 
itive constant. The Schmidt gap AA can be expressed 
as: 



AA(£) = Ai - A2 = 



1 - 
ZR{q) 



1 



lc/12 



(3) 



where we exploited the fact that the largest eigenvalue 
of the reduced density matrix is directly related to the 
single copy entanglement 5*1 = — In Ai [33| and goes as 
Ai ^ £~'^/^'^. The coefficients afc > have the form 
A + iV {N > 0) where {A} are the scaling dimensions of 
the primary operators of the conformal theory describing 
the critical point and {N} denotes the set of integer pos- 
itive numbers labeling the conformal levels. Each level 
has a computable degeneracy d!]^\ The numbers A + N, 
defining an ordering labeled by the index k above, cor- 
respond to the eigenvalues of the operator Lq [1^, HI] . 
Since the eigenvalues of this operator arc fixed by the 
central charge so is the Schmidt gap at criticality for a 
finite size chain. 

Out of criticality, but still in its vicinity, the entangle- 
ment spectrum is expected to acquire specific features, 
loosing the universal dependence on the central charge 
valid at the critical point. It can be inferred that for an 
integrable model, the trace of the nth-power of reduced 
density matrix pn for the semichain described close to 
criticali ty j 2]l| can be expressed by a character sum as 

MMm-. 



Ea «A(g) XA(g") 



where XA (q) are the conformal characters |26l . 134 



q 



(4) 



(5) 



N=0 



Indeed, in these models the eigenvalues of are all of the 
form q°'^+^ and with computable degeneracies d^'' . The 
parameter 5, depending now on the correlation length ^, 
tends to 1 at criticality — >■ 00), there it agrees with the 



modular parameter q = exp[— 27rK/ log(^/77)], appearing 
in ([21) , in the limit £ — > 00 . 

This observation allows us to express the Schmidt co- 
efficients near criticality as: 

>^k = ^ — — = , (6) 



EaOaI?) XA{q) 



ZR{q) 



Zji[q) is added as a trace normalization for the {Afc}. 
The difficulty in evaluating the above expression comes 
from the lack of knowledge of the coefficients oaC?) ap- 
pearing in (131) whose values are only known at criticality 
where they fulfill a a (9 = 1) = 1- This limit can be 
inferred from the fact that Trp^ is expressible at crit- 
icality just simple sum of the characters, with all 
coefficients equal to 1 [13: 113 • For continuity the a/^{q), 
when nonzero, should be close to 1 in the scaling limit. 
For sake of generality we assume that they all have a q- 
dependence. Actually, the fact that on a certain side of a 
given phase transition some a{q) vanish for every q seems 
to suggest that each a[q) generally assumes only values 
or I, probably due to some protecting symmetries, the 
same ones underlying integrability [1^ . We leave this is- 
sue as an open point. Finally, notice that the positivity 
of Afc implies positiveness also of a^(q). 



A. Transverse magnetic field Ising Model 

We will go through the detailed calculation for a known 
case, the spin-1/2 transverse Ising model, having central 
charge c = ^ at the critical point. Obviously, the same 
discussion holds for all the lattice models that at critical- 
ity belong to the same universality class. A direct anal- 
ysis of the entanglement spectrum obtained numerically 
(see subsection II Cp shows that in the low-temperature 
(ferromagnetic) phase the character entering in ^ is 
Xi('z) (indeed the entanglement spectrum is doubly de- 
generate) while in the high-temperature (paramagnetic) 
phase Xo{q) and X^iq) are relevant. We consider the sec- 
ond case, since it also reproduces the Ising transition in 
the spin-1 model approaching the Haldane phase from 
the Ncel one. Setting n = 1 in (jlj, the Schmidt gap 
AA = Ai — A2 can be written as: 



AA = 



q 



aoiq) - ai{q)q2 



l-/3o.i('7)g = 



ZR{q) 



(7) 

where /3g i((?) — ai/2{q)/aQ{q) and we have that Si = 
— InAi ss In^. In this particular model 00(9) and a 1 (q) 
can be probably computed expressing pn as a thermal 
ensemble of free fermions [36j , however we will see that 
the precise expressions for them are not strictly required 
in our calculation. 

We want to evaluate q in terms of ^ in the limit q — > 1, 
i.e. close to the critical point, where unfortunately the 
series in Eq. (0) it is not so useful since it must be re- 
summcd before taking the limit q — >■ 1. However, the 



4 



Virasoro characters transform hnearly under a modular 
transformation q q: 



XA{q) 



(8) 



A' 



where q and q 
relation |25|. ul 



are linked each others by the standard 



Inq \nq = 47r 



(9) 




The scries ^ expressed in q has the virtue that the criti- 
cal limit g — > can be taken before to resum all the terms 
[25| . is called modular matrix and characterizes the 
modular transformations on the torus for a given CFT 
[1^. For the Ising CFT (c = i) it reads 



(10) 



in the basis (xo(<z)j X^{<l)^ X-i^il))'^ or equivalently 
(Xo(g), X^il), X-ig{Q)V (the matrix is nilpotent of or- 
der two). 

Since the first column, concerning the expression of the 
Xiil) in terms of Xo(9)i contains only positive numbers, 
as well as the coefficients a^iq) in ([4]), it is clear that, af- 
ter the transformation ([5]), the expression in q for Trp^, 
analogous to (|4]), surely contains the character Xo{<j)- No- 
tice that this claim does not depend strictly on the par- 
ticular expression (j4|). We consider now the functional 
dependence of Tr in terms of q. It is straightforward 
to extract the leading behavior in the limit g — >■ 0, as well 
as the first corrections. The leading term comes from the 
character xo(9) (surely present in Trp^(g)): 



ao{q)S^o+ai{q)Sl 



(11) 



which can be used to obtain the Renyi entropy 
1 



1-n 



(12) 



-^h+i)lng + ln ao(g)^S + «i (9)^4 



Comparing (jl2p with the leading term for the Renyi en- 
tropy out of criticality [3 [25| 



Sn 



C 

12 



InC, 



(13) 



one arrives to q 
q = exp 



^ {v being a constant) and 



lnr)-21n4 

Finally, the term 



In 



is related 



a„{q)S^ + ai{q)Sl 

to the Afheck-Ludwig boundary entropy g [37[. Notably 
a measure of it can give direct information on ao(g) and 
ai (g). Since aA(l) = 1, g^ is expected to be well approx- 
imated by g"^ « In (^Sq 



5? 



Inserting the expression for g in ([71) and expanding the 
exponential yields: 



AA(e) 



5T5(lni)-2 10 



^T2 



CTTI(ln))-210 

(14) 

where c(g) = 1 — /Jq |(g)- The limit c(^) — > is not 
trivial and one cannot naively approximate c(^) w 0, the 
peculiar functional form of c(^) determines instead which 
term in the r.h.s of is dominant for ^ — >■ 00 (critical- 
ity). 

Symmetry arguments discussed at the end of the last sub- 
section and general arguments from the standard scaling 
theory (we take into account that no marginal operators 
are present) [s^ lead us to believe that c(^) — )• more 
rapidly than . Moreover, if our conjecture that aA{q) 
are always either or 1 is true, this term is exactly van- 
ishing. Using the scaling relation £, = \g — 5c|~^ and 
c = ^, we arrive to the final expression: 



AX{\g-g,\) = A 



\9- 9c\ 



B - ln|.g 



(15) 



with A = 7r^(47r)24 « 10.97. The value for B instead is 
not simply derivable. The factor required for the es- 
timation of A, can be derived from the exponential decay 
rate of the two-point correlation function of the energy 
density operator [2l| and considering that the mass gap 
scales as: M = 2TT\g — gc\ [1^. Notice that in p5|) . the 
dependence on the precise characters content is only in- 
cluded in the parameters A and B. The relation (fT4)) . in 
particular, is expected to be universal (but with varying 
c(g) and 69 1(9)) for integrable models. 
Formula (fT?|) suggests that the scaling of the Schmidt 
gap depends mainly (in its functional form) on the cen- 
tral charge and on the correlation length of the system 
(in turn depending on the Hamiltonian operator per- 
turbing the critical point), ingredients already encoded 
in (fT3|) . The precise expressions for A and B depend 
instead on the Schmidt eigenvalues, via the parameters 
a\{q). However the derivation presented in this section 
works only for integrable perturbations and we leave as 
an open question whether in a general non-intcgrable case 
a stronger dependence on the full entanglement spectrum 
holds. 

Notice that for g gc the expression (fT5|) for the 
Schmidt gap tends to: 



AA(|ff-5c|) = -^ 



— 24 



ln|5 - 9c 



(16) 



to be considered the correct leading term in the thermo- 
dynamic limit. The factor ^, \_ — r correcting the power 

law was already noticed in [2 a . |40j . 

However we point out that, dealing with finite chains, 
the B term in the denominator of Eq. (|15p is required to 
be kept, since the quantity In^ in (fTH) is prevented by 
the finite size to get sufficiently big in magnitude such 
to make B negligible. For this reason the comparison 
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with DMRG data in the next section wiU be performed 
assuming (|15p . 



B. Longitudinal magnetic field Ising model 

In this case the character content in ^ can be guessed 
by simple considerations. In a generic case, the expres- 
sion (m at the critical point contains all the characters 
of the CFT describing it, while outside the critical point 
only some of them appear on each side of the transition. 
It's then clear that the union of the sets of characters ap- 
pearing in Trp^(q) on every side of the transition must 
be equal to the entire set of characters appearing in (j4|) at 
criticality. The Ising model with longitudinal magnetic 
field displays a QPT between two phases with no residual 
symmetry, thus no constraints on the character content, 
directly related to the degeneracies of the Schmidt eigen- 
values. The two phases have both a single vacuum, they 
share the same S'-matrix and they can be mapped into 
each other 2l[. For this reason we conjecture that, out 



of criticality, Trp^(g) contains all three characters: 



^ «o(g)xo(g") + aj,(g)x^(g") + a.(g)xi(g") 

(17) 

where again a^{q) > 0. Notice however that in order 
to obtain a Schmidt scaling law analogous to (fT5|) . the 
explicit character content is not explicitly required, again 
thanks to the positivity of the coefficients aAil) ^^nd to 
the peculiarity of the matrix ()10p . Setting again n = 1 
in (HI and ([TT]) and taking into account (O, the Schmidt 
gap can be expressed as: 



AA 



ao{q) 



ZR{q) 



(18) 



Using the same arguments and calculations as in the pre- 
vious subsection and the scaling relation ^ = a \g — gd^"^ , 
we obtain the formula: 



AA(|g-5c|)=A^ 



B - ln|.g - 5c 



(19) 



with A = 19.27 (while again the value for B is not easily 
derivable) and in the thermodynamic limit 



AA(|.g-5c|) = -i"^ ^' 



— (7„ 45 



ln|5 - 9c 



(20) 



We explicitly inserted in (|T9|l the values of the sealing 
exponents ly ~ 8/15 and a ss 0.38, this last values was 
estimated numerically in [4l| from the decay rate of the 
two-point correlation function for the spin operator. 
In again the dependence on the precise characters 
content in ([TT]) is only in the parameters A and B. 



C. Comparison with numerical data 

At this point we want to check the predictions out- 
side criticality obtained analytically against the results 
obtained by the scaling of the Schmidt gap when solving 
exactly or numerically the corresponding spin model: 



H 



Isirig 



The model ([?!]) reduces (i) for Bj. ~ to the transverse 
Ising model, that is critical at J = B^, and (ii) to the 
longitudinal Ising model for B^/B^ ^ 0; both models 
are described at criticality by a CFT with central charge 
c = 1/2, but they are in different universality classes 
j2ll [26j . The entanglement spectrum and the Schmidt 
gap for both models were investigated in [l3|. For the 
transverse Ising model the entanglement spectrum can 
be computed employing the Jordan- Wigner and Bogoli- 
ubov transformations to map the model into a system of 
non interacting fermions, so that the Schmidt eigenvalues 
can be straightforwardly extracted [11] . The Ising model 
with a longitudinal field perturbation is not integrable in 
the lattice and therefore has to be solved numerically, by 
means, e.g. of DMRG at J = and for different values 
of Bx/Bz 7^ 0. This model is however integrable in its 
infrared continuous limit, as well as the transverse Ising 
chain. In both cases we found [l3] that, when approach- 
ing the critical point, the Schmidt gap displays scaling 
behavior with the size of the system i and it allows for 
a FSS analysis, as for a standard local order parameters 
in a second order phase transition (28j : 



AA(£,g)^r^-/VAA(|5-5cK'/'') 



(22) 



The quantity v characterizes the divergence of the cor- 
relation length, while /3aa plays the role of an "order 
parameter" exponent. Indeed, the ansatz (j22p implies a 
power-law hypothesis for the Schmidt gap in the infinite 
limit: 



AA 



(23) 



While for local order parameters as the niagnetization 
this hypothesis is under control (see note [42]), for the 
Schmidt gap, being a non local quantity, the justification 
is highly non trivial, and it must be checked a posteriori, 
for instance probing the correctness of the v exponent 
obtained from (|22|). 

From the scaling of Schmidt gap by p2|) we extract the 
following results (see [13): for the transverse Ising model 
the critical point is located at J/S^ = 1 and the critical 
exponents /3aa = 0.124±0.002 and l.OOiO.Ol. These 
values agree very closely with the critical exponents on 
the Ising transition, suggesting that the Schmidt gap in 
the vicinity of the critical point correctly signals the lo- 
cation of it and it scales universally with critical expo- 
nents related to the CFT describing the transition. For 
the longitudinal Ising model the scaling of the Schmidt 
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gap leads to the critical exponents: v = 0.50 ± 0.05 and 
/3aa = 0.055 ± 0.005, which again are in good agreement 
with the exact values of the corresponding exponents for 
the mass gap (y = 8/15) and for the magnetization order 
parameter (1/15). 

To test the versatility of our estimates we compare them 
with the numerical results for the Schmidt gap in the 
transverse and longitudinal Ising model with fixed system 
size. Surprisingly enough, the scaling behavior obtained 
from the power law ansatz and from the expressions 
([T5|l and ([T9|) are strongly compatible, as displayed in 
Figs. O and (HI). 

For the transverse Ising, the fitted parameters in (|23p are 
/3a A = 0.114 and k = 1.05. Assuming instead (fTSjl. the 
fit yields B = 8.94 and A = 10.27, to be compared with 
the expected value A = 10.97. The succesful comparison 
between numerical data and the derived scaling law (flSl) 
represents a further ar gum ent in favor of the Baxter and 
Cardy's conjecture [20n25l. [26| . Moreover it supports the 
claim that for finite size systems the B term in (jl5p is 
needed. 

For the longitudinal Ising model, the power-law fit ([^5]) 
of the numerical data for a chain of 1536 sites gives 
/3aa — 0.03 and k = 1.16; assuming instead (fT9|) we 
obtain S = 71 and A ~ 84.2, to be compared with a 
theoretically expected value A = 19.27. Compared with 
the transverse Ising model, the larger error between the 
measured and the expected values of A is probably due 
to the smaller size of the analyzed chain and to the bigger 
sensitivity of it to finite size corrections. 
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FIG. 1: Schmidt gap closing in the vicinity of the critical point 
for the transverse Ising model. Comparison between the ana- 
lytical prediction from CFT (|15|) (solid grey line), the power 
law scaling (|23|l (bin dashed line) and the numerical values 
(crosses) obtained by exact diagonalization. We consider a 
chain of 12000 sites. 
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FIG. 2: Schmidt gap closing in the vicinity of the critical 
point for the longitudinal Ising model. Comparison between 
the analytical prediction from CFT (|19|l (solid grey lines) , the 
power law scaling (|23|l (bin dashed line), and the numerical 
values (crosses) obtained by solving the model with DMRG. 
We consider a chain of 1536 sites. Notice that the fit is re- 
stricted to the six largest values. 



spoiled here by the infrared relevant perturbations (in the 
renormalization group (RG) sense) of the critical theory. 
One could expect that, even for non-integrable cases, very 
close to the critical point {q ~ 1) or to an integrable tran- 
sition this effect to be very small and the conjecture to be 
still approximately valid. However this regime is proba- 
bly too small to allow us a reliable study of the scaling of 
the Schmidt gap and the numerical analysis is the unique 
way to proceed. 

In what follows we analyze the scaling properties of the 
entanglement spectrum in the various quantum phases 
appearing for the spin-1 chain in the presence of an uni- 
axial anisotropy D. Some of the cases here discussed are 
non-integrable to the best of our knowledge. The Hamil- 




e/Tt 



II. SCHMIDT GAP IN SPIN-1 CHAINS: 
SCALING IN ABSENCE OF INTEGRABILITY 

For non-integrable models the conjecture used to de- 
rive the scaling of the Schmidt gap in the previous sec- 
tion does not hold, since the conformal degeneracies are 



FIG. 3: Phase diagram of the spin-1 Hamiltonian H24|) as a 
function of 9 and D computed by DMRG. The details about 
the phases can be found in 30]. We also put an arrow and 
a corresponding number for the four critical points that we 
study more extensively in the text. 
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tonian describing a spin-1 chain is given by 

H^Y.^m+^Y.^^-^^ (24) 

i i 

where Si = {Sxi, Syi, Szi) are the i-th site angular 
momentum operators and Hi{9) ~ cos{d) Si ■ Si-^-i + 
sm{6){Si- Si+i)'^ . For D = 0, the model is the well known 
bilinear-biquadratic spin-1 chain. The phase diagram for 
this chain as a function of G [— tt; tt] is known, see for 
example [i^. For — 7r/4 < 9 < 7r/4, the system is in the 
Haldane phase. At 9 = — 7r/4 we obtain the so called 
Takhtajan-Babujan spin-1 chain, a model integrable on 
the lattice. From the solution of the Bethe Ansatz equa- 
tions, the only elementary excitations are known to be a 
doublet of gapless spin-1/2 spinous, with total spin or 
1 nil. A detailed description of the phase diagram for 
D ^ and its realization by ultra-cold gases is reported 
in |29l . [soj . As shown in Fig. [31 the uniaxial anisotropy 
leads to several phases surrounding the Haldane one. 
Those are: (i) the Neel phase, appearing for negative 
D; (ii) the dimer phase, present for negative D, but also 
for very small positive values of D; (iii) a large- 1? phase 
whose ground state is connected adiabatically to the state 
in which all spins have a zero angular momentum in the z- 
component; (iv) a critical gapless phase. Further to these 
phases, the model exhibits a ferromagnetic phase, an XY- 
ferromagnetic phase and a critical XY-antiferromagnetic. 
To characterize the Haldane topological phase and to lo- 
cate precisely its boundaries, we compute the dimer order 
parameter V ~ Hi{9) — Hi^i{9), the staggered magneti- 
zation per site = 1/i "Yliii^^Y Szi, acting as an order 
parameter for the Neel phase, and finally the entangle- 
ment spectrum. Additionally, the Haldane phase is char- 
acterized by a non vanishing string order parameter (43| : 

0= lim (5^,; exp[i7r SzASzi+r) ■ (25) 

j='i+l 

A. Neel-Haldane phase transition 

As a first step we concentrate our analysis on the 
Neel-Haldane phase transition along 6 — the uniaxial 
anisotropy D is now the control parameter, as schemati- 
cally shown in Fig. [3] with arrow 1. In our DMRG sim- 
ulations we apply a small magnetic field to the first spin 
to select one of the two degenerate ground states. The 
FSS analysis for the staggered magnetization yields val- 
ues for the location of the critical point and for the crit- 
ical exponents respectively Dc = —0.315, /3 = 0.11 and 
V — 1.01, very close to those obtained from the FSS of 
the Schmidt gap Dc = -0.315, P = 0.11, v = 1.04. These 
results arc comparable with the ones from the quantum 
MonteCarlo ^ = -0.316, /3 = 0.147, i/ = 1.01 and 
are in agreement with the conjecture [45| that this tran- 
sition belongs to the transverse Ising universality class 
/3 = 0.125, V — \. The 12% discrepancy in our esti- 
mate of /3 can be attributed to the relatively small size of 



the considered chains, the same for MonteCarlo results. 
This similarity suggests us to compare the statistics of the 
Schmidt eigenvalues from DMRG with that one obtained 
by the character analysis performed out of criticality for 
the transverse Ising model in Sec. IH By direct inspection 
we checked that the statistics in the two cases perfectly 
agree, confirming that these two different models are in 
the same universality class and therefore share the same 
entanglement spectrum outside, although still close to, 
criticality. 

As for the transverse Ising model, the scaling of the 
Schmidt gap leads to similar critical parameters as those 
found from the FSS of the staggered magnetization. This 
correspondence seems to be closely related to the obser- 
vation that, even far from the transition point, staggered 
magnetization and Schmidt gap behave almost identi- 
cally, as shown in Fig. [H 
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FIG. 4: (Figure from j3Cl]) The Schmidt gap, the string order 
parameter and the staggered magnetization for the Hamilto- 
nian H24|) on the 6 = line, shown as a function of D/ J . The 
shaded area represents the scaling regions close to the critical 
points where our numerical analysis is less reliable. 



B. Dimer-Haldane phase transition 

We analyze now the transition from Dimer to Haldane 
phase along D ~ Q and changing 9 (see arrow 2 in Fig.[3l). 
The conformal point related to this transition is the pre- 
viously mentioned Takhtajan-Babujan point [43 | and it is 
known to be described by a SU (2) current algebra (Wess- 
Zumino-Witten, WZW) theory at level 2 (following the 
standard notation we will denote it as SU {2)2)- This the- 
ory is characterized by a central charge c = 3/2 and two 
primary relevant operators with spin 1/2 and 1 respec- 
tively and conformal dimension 3/16 and 1/2 [21]. 126} . 
We focus on the FSS of the dimer order parameter V 
and the Schmidt gap AA. The lattice model ([24|l with 
D ~ Q \s SU{2) invariant, therefore, provided that the 
cut in space does not break explicitly this symmetry jsH , 
we expect the Schmidt spectrum to organize in multiplets 
corresponding to the representations of this group. For 
every considered even size L, the cut in the middle point 
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breaks an S'?7(2)-singlet either in the dimerized or in the 
Haldane phase, depending of the parity of L/2. How- 
ever the effects of the SU{2) breaking are negligible for 
long chains, as it usually happens for boundary effects, 
then at the end we expect the spectrum to be always ar- 
ranged in multiplets, no matter the position of the cut. 
This picture is confirmed in Fig. [51 where multiplets 
are clearly visible. We notice also the peculiar distribu- 
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FIG. 5: Schmidt spectrum close to the Dimer-Haldane phase 
transition for D = as a function of 9 and for L — 384. 
The critical point is at 9/n = —1/4 and it is indicated by a 
vertical dashed line. The discrepancy between the location 
of the critical point and the zone where the statistics of the 
entanglement spectrum changes is a finite size effect. We 
indicate the degeneracies with small clusters of dots in the 
two phases. Lines joining the numerical data are only a guide 
to the eye. 



tion of the Schmidt eigenvalues in both the Haldane and 
the dimer phases: the first one has only even multiplets, 
while the second one contains only odd multiplets. This 
structure agrees with recent results obtained in [i^, \^ , 
where all the possible phases of the ID gapped and short- 
range interacting Hamiltonians are classified by assum- 
ing a matrix product state (MPS) ansatz for them and 
studying the isometries of MPS. There a correspondence 
was found with the different projective representations 
of the local invariance group G for the analyzed energy 
eigenstate (generally the ground state) . It is then not dif- 
ficult to see that one can equivalently classify the linear 
representations of the covering group G in subsets that 
are one-to-one related to its center group (the precise re- 
lation comes directly from group representation theory, 
see di]). The Hamiltonian ([H]) along the line D = 
is invariant under SO{3) = SU{2)/Z2 therefore different 
phases corresponds to the two sets of linear representa- 
tions of SU{2) labeled by the different elements of its 
center Z2'. the spin integer and semi- integer represen- 
tations. This link translates directly in the content of 
multiplets of the Schmidt spectrum, as shown clearly in 
Fig.El 



The FSS analysis for the dimer order parameter V and 
the Schmidt gap AA is shown in Fig. [6l The results 
for the critical point extracted from 2? and AA are 
9c = — 0.2497r and 9c = — 0.2547r respectively, very close 
to the exact result 9 = — 7r/4, confirming the effective- 
ness of the Schmidt gap as an order parameter signaling a 
phase transition. The estimates for the correlation length 
exponent arc v — 1.41 and v = 1.45 respectively and they 
agree with the exponent calculated directly from the scal- 
ing of the energy gap: i/ = 1.43. Again the scaling law 
ansatz ([22]) turns out to be effective, yielding the correct 
value for v. Finally, the critical exponent for the dimer 
order parameter is /3 = 0.66 and for the Schmidt gap is 
P = 0.20. 

Let us try to gain some insight on these exponents us- 
ing arguments from CFT. The D = line has SU{2) 
symmetry, therefore the most general perturbation of 
the critical point SU{2)2 along this line is expected to 
be a combination of certain powers of the singlets con- 
structed by the relevant fields of the SU{2)2 WZW the- 
ory. The singlet constructed by the spin-i multiplet is 
{cl>i4>^i ~ 4'^i'Pi)j while the singlet constructed by the 

2 2 2 2 

spin-1 muhiplet is (0^0o - 0f 0-i - (f>f4i^i)- Notice 
that here R and L denote the chiralities of the confor- 
mal SU {2)2 algebra [H Hlj. The last possible singlet 
is the scalar product of the chiral and anti-chiral SU (2) 
vectorial current sj = Jr ■ Jl- Eventually we obtain 
that the most general relevant perturbation along D = 
i*^ Z]t, (*i/2)" + (si)'" + Y.ki^j)''- Each (n, m, fc)-powcrs 
can be further decomposed in a sum of terms derived 
from the one point expansions (OPE) of the (n, m, k)- 
product of chiral field s I26l. An explicit calculation em- 
ploying bosonization j26l l49l [soj shows that in our case 
the Hamiltonian density of the continuous theory close 
to the critical point is 

H = Hsu{2h+a{9)si+bsj, (26) 

where a{9) oc 9 — j. Interestingly, si turns out to be 
a mass term for the three Majorana fermions describing 
the SU(2)2 point. The application of both si and sj 
was argued to give rise to a path from SU{2)2 to the 
Heisenberg point at 6* = 0. described by a sigma model 
without topological term [40|. Since both a relevant op- 
erator and a marginal one appear in ()26p . then a simple 
power law behavior for the mass gap and the magneti- 
zation is not expected even in the infinite limit, but a 
logarithmic correction to it. This picture was confirmed 
in j5l| . Closely related, for a finite chain a marginal op- 
erator has a dramatic effect on the scaling dimensions of 
the operators of the infrared continuos theory, as argued 
in [SOj. In particular, one can define effective scaling di- 
mensions which fiow logarithmically slowly in L to their 
asymptotic values. For the dimer order parameter and 
for a chain havin g 25 6 sites a scaling dimension x = 0.48 
was measured in [50| , to be compared with the infinite L 
value a; = 3/8 = 0.375. The value 0.48 is very close to 
the value x = 0.47 that we obtained for a chain of the 
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FIG. 6: Finite size scaling close to the Haldane-dimer transition. From the left size, the first and second panels are the FSS 
results for the Schmidt gap AA, while the third and fourth panels refer to the FSS of the dimer order parameter. 



same length L = 256 using the relation x — ^ and the 
measured values (3 = 0.66 and = 1.41. 



point 6 ~ 
= 0.75. 



-0.457r and to the exponents ly = 1.0 and 



C. Large D - Haldane phase transition 



III. DISCUSSION AND CONCLUSIONS 



We consider now the transition between the large D 
phase and the Haldane phase. This phase transition line 
is described by a CFT with central charge c = 1 related 
to a free boson compactified on a circle (H^]. The bo- 
son velocity, the radius of the circle and importantly the 
operator content of the theory depend on the parame- 
ters 9 or D. This fact implies that the critical exponents 
are continuously varying changing the transition point 
along the critical line [52[. We focus on the transition 
at 6* = (sec arrow 3 in Fig. [3]). Our results for the 
Schmidt gap provide the critical point at D^^ = 0.96 
and u'^'^ = 1.56. These results are again compatible from 
those obtained with Quantum MonteCarlo [1^ simula- 
tions: D9^^ = 0.971 and v'^^'^ = 1.4. Also in this case 



the scaling law ansatz (|22|) looks to be pretty accurate 
and provides a critical exponent /3aa = 0-3. Again the 
Schmidt gap turns out to be a reliable order parameter. 



D. Neel - dimer phase transition 

For completeness we investigate now the QPT between 
the dimer and the Neel phases (arrow 4 in Fig. [3]) and an- 
alyze the scaling of the Schmidt gap along this transition. 
Our DMRG calculations of the entanglement entropy S 
show that, decreasing the parameter D from D = to 
Z? = — oo, the central charge along the transition line 
varies continuously from c = | to c = 1. In order to ex- 
tract c we used the formula S" = | Log L describing the 
sca ling of the half chain entanglement entropy at critical- 
ity [12, [3 • The observed peculiar behavior it is possible 
because the (1-f 1)D conformal field theories have a con- 
tinues set of central charges and Verma modules for c > 1 
[2i,[3l. 

Also in this case we expect the Schmidt gap to close at 
the critical point. Indeed, for D = —0.9, AA exhibits 
critical scaling around 6 ~ — 0.457r with critical expo- 
nents (3^^ = 0.19 and v = 1.2 while the FSS analysis 
of the dimer order parameter leads to the same critical 



Here we have analyzed different QPTs, either intc- 
grable or not, by the scaling properties of the Schmidt 
gap, i.e. the difference between the two largest non 
degenerate Schmidt eigenvalues We have derived 

analytically the scaling of the Schmidt gap for the 
spin-i Ising chain in a transverse or a longitudinal 
magnetic field and we have checked the effectiveness of 
these formulae on the numerical data obtained by direct 
diagonalization or DMRG respectively. This succesful 
comparison, especially in the case of the transverse Ising 
model, supports the old conjecture by Baxter [23| and 
Cardy [2g relating the Schmidt eigenvalues close to 
a phase transition to the characters of the conformal 
theory at the critical point. 

By means of DMRG simulations, we have studied the 
Schmidt gap for various phase transitions of the bilinear- 
biquadratic spin-1 chain with a quadratic Zceman term. 
This study allows us to extend some of our results to 
a wider set of QPTs than the transitions of the two 
Ising universality classes, described for instance by 
non-integrable infrared effective models or characterized 
by the presence of continues symmetries. 
In particular, our findings have shown that the Schmidt 
gap correctly signals a QPT and, approaching the critical 
point, it displays a scaling behavior as a conventional 
local order parameter. This scaling is characterized by 
critical exponents closely related to the CFT describing 
the critical point and to its specific perturbations. 
Moreover, thanks to the power of the FSS technique, 
the mass scaling exponent v turns out to be computable 
with good precision from the study of the Schmidt gap. 
This allows us to identify the Schmidt gap as a non-local 
order parameter, particularly valuable whenever a 
"standard" local order parameter is not known, for 
instance because it cannot be defined at all, as in the 
case of topological phases, or because the nature of the 
phase itself is obscure. 

At first glance, the effectiveness of the Schmidt gap to 
locate phase transitions would stem from the behavior 
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of the entanglement spectrum, that at criticality and 
in the thermodynamical limit must tend to a continuos 
distribution: this simply implies the closure of the 
Schmidt gap. Although finite size effects generally 
prevent this closure, the scaling behavior survives and it 
can be recovered by FSS analysis. 

The analytic derivation of the Schmidt gap scaling for 
transitions in the Ising universality classes suggests that 
a scaling behavior should also hold for the difference 
between higher order Schmidt eigenvalues |Ai — Aj|, how- 
ever in the numerical data from exact diagonalization 
(for the transverse Ising model) and form DMRG no 
scaling and not even closing is observed for them. We 
leave as an open point whether this lack of scaling is 
only due to finite size effects. 

Remarkably, one result of the present paper is also a 
direct check of the results in [ia, \^ about the classi- 
fication of topological order in short-range interacting 
gapped quantum many body systems: the analysis of 
the entanglement spectrum along the parameter line 
D = 0, including the vicinity of the dimer-Haldane phase 
transition, displays a multiplet structure with statistics 
in agreement with these results. 

The observed properties of the entanglement spectrum 
and of the Schmidt gap can be an important boost for 
future developments concerning the characterization of 
strongly interacting systems with higher dimensionality/ 



and or with long range and order, as in ultracold polar 
molecules systems [6| or in presence of deconfined 
quantum critical points 0, 01 , where moreover a peculiar 
entanglement structure has been recently found (53j . 
Indeed, first results showing scaling properties of the 
entanglement spectrum near QPT for 2D quantum sys- 
tems have very recently appeared [3 13 • The analysis 
we propose is particularly relevant for fermionic systems, 
when the sign problem forbids trustable numerical results 
obtained via a MonteCarlo approach. Some preliminary 
attempts to describe long-range interacting systems in 
one dimension are already available [53 |. A central topic 
along this direction is the classification of the topological 
order and phases in such systems, enlarging the logical 
scheme adopted in [i^, |43]. The entanglement spec- 
trum has been shown to encode important features of 
topological order and its analysis, along the lines of the 
present paper, is expected to have a primary importance. 
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